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We use the Glauber- Gribov multiple scattering formalism and the theory of leading twist nuclear 
shadowing to develop a method for the calculation of leading twist hard coherent diffraction in 
' hadron-nucleus processes. We demonstrate that soft multiple rescatterings lead to the factorization 

o ' 

breaking of hard diffraction in proton-nucleus scattering, which is larger than in hadron-nucleon 
scattering. At the LHC and RHIC kinematics, we compare the hard diffractive to e.m. mechanisms 

in 
o 



of hard coherent production of two jets in proton-nucleus scattering. We study the xjp, /3 and A- 



' dependence of the ratio of the dijet production cross sections due to the two effects, R. We 

a: 

demonstrate that in proton-heavy nucleus hard coherent diffraction at the LHC, R is small, which 
offers a clean method to study hard photon-proton scattering at the energies exceeding the HERA 



energies by the factor of ten. On the other hand, the use of lighter nuclei, such as Ca, will 
allow to study the screened nuclear diffractive parton distribution. Moreover, a comparison of the 
dijet diffractive production to the heavy-quark-jet diffractive production will estimate the screened 
nuclear diffractive gluon PDF, which will be measured in nucleus-nucleus ultraperipheral collisions 
at the LHC. 
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I. INTRODUCTION 



In hadron-hadron scattering at high energies, diffractive processes are characterized by 
the rapid t-dependence and by the absence of detected particles in a certain region of the 
final phase space or, in other words, by the presence of the rapidity gap. When a hard scale is 
present in diffractive scattering, such processes are called hard diffractive. The phenomenon 
of hard diffraction in proton- antiproton scattering was first discovered in the pp — > pX 
reaction at the the SPC collider at CERN, when it was observed that the diffractive final 
state X with the invariant mass in the range 105 to 190 GeV contained jets with the 
transverse energy between 5 and 13 GeV (the hard scale is given by the jet transverse 
momentum) [lj. Later, hard diffraction in proton-antiproton scattering was studied at the 
Tevatron in dijet, W, 6-quark and J/ty production, see ^ for a review. The cross section 
of each hard diffractive channel constitutes approximately 1% of the contribution of the 
corresponding channel to the inclusive pp cross section. 

Turning to electron-proton deep inelastic scattering (DIS), it was one of HERA suprises 
to observe that hard diffractive events characterized by a large rapidity gap between the 
virtual photon and the proton fragmentation regions constitute about 10% of the total rate 
of events [3]. In DIS, the hard scale is provided by the virtuality of the photon, Q 2 . 

In the theoretical treatment of soft diffractive processes, the key role is played by the 
concept of the "Pomeron", the Regge trajectory with vacuum quantum numbers, which 
provides the diffractive exchange and determines the high-energy behavior of elastic and 
diffractive scattering amplitudes. In the context of hard diffraction, the notion of Pomeron 
appears as follows. The QCD factorization theorem for hard diffraction in DIS |4| enables 
one to introduce universal diffractive parton distribution functions (PDFs), which can relate 
such processes as inclusive diffraction, dijet diffractive production, D*-meson diffractive 
production, etc. Making an assumption that the diffractive PDFs can be factorized into the 
product of two terms, representing the Pomeron flux and the Pomeron parton distributions, 
one can effectively study the parton content of the Pomeron, similarly to the parton content 

fin nnri 

of the nucleon in inclusive DIS [5|, |6|, |7|, |8|, |9j. 

A comparison of hard diffraction in proton-antiproton scattering at the Tevatron to hard 
diffraction in electron-proton scattering at HERA indicates the breakdown of the QCD fac- 
torization: The use of diffractive PDFs extracted from the HERA measurements significantly 
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101 ]. This can be explained 



overestimates the rates of hard diffraction at the Tevatron {2], 
by the absorptive effects associated with multi-Pomeron exchanges, which make the gap 
survival very unlikely in the case of hadronic collisions 11|, or by the gradual onset of the 
so-called black disk regime for the proton-proton collisions at small impact parameters 12]. 

Turning to diffraction in hadron-nucleus and lepton-nucleus scattering, the situation can 
be briefly summarized as follows. Soft coherent (without the nuclear break-up) diffrac- 
tion in hadron-nucleus scattering at high energies can be successfully described within the 
framework of the Glauber-Gribov approach by taking into account cross section (color) 



n 



fluctuations in the hadronic projectile [l3|, |l4|, Ilia, Il6|, see also Sect. HI1 



n 



In DIS on nuclear targets, nuclear diffractive PDFs at small values of Bjorken x can be 
expressed in terms of the nucleoli diffractive PDFs. which arc known from the HERA stud- 
ies 17|. This approach to nuclear diffractive PDFs and to usual nuclear parton distribution 



functions is based on the Gribov's connection between the nuclear shadowing correction and 
the elementary diffractive cross section, the QCD factorization theorem for hard diffraction 
in DIS |4j and the QCD analysis of HERA data on hard diffraction [JflfllslB, see Q 
and references therein and also Sect. [Ill 

In this work, we study hard diffraction in proton-nucleus collisions. As an example, we 
derive the expression for the cross section of the hard coherent diffractive production of two 
jets in proton-nucleus scattering. We analyze the suppression of hard diffraction in proton- 
nucleus scattering compared to hard diffraction in proton-proton scattering. A numerical 
analysis of the corresponding suppression factor enables us to quantify the QCD factorization 
breaking in hard proton-nucleus diffraction due to the soft screening (absorption). 

The paper consists of two parts. In Sect. [Til we derive a general expression for the 
nuclear modifications of hard diffraction in proton-nucleus scattering. We demonstrate that 
because of multiple soft rescatterings, hard diffractive processes such as production of two 
jets, heavy flavors, etc. are suppressed at the LHC and RHIC energies stronger than soft 



inelastic diffraction, which in turn is expected to be strongly suppressed [13|, [l5|, [16] 



In Sect. IIII[ we compare the contribution of hard coherent proton-nucleus diffraction into 
two jets (including heavy-quark jets) to the e.m. contribution, when the final state containing 
two hard jets is produced by the coherent nuclear Coulomb field. We demonstrate that the 
e.m. contribution dominates proton-heavy nucleus (such as 208 Pb) scattering at the LHC, 
which provides essentially a background- free method to study very high-energy jp scattering 
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at the LHC through ultraperipheral proton-nucleus scattering. We show that using lighter 
nuclei, which do not produce such a strong flux of equivalent photons as 208 Pb, one can 
study screened nuclear diffractive PDFs. In this case, a comparison of the dijet diffractive 
production to the heavy-quark-jet diffractive production will measure the nuclear screened 
diffractive gluon PDF. The latter can be compared to the nuclear diffractive PDFs, which 
will be measured in nucleus-nucleus ultraperipheral collisions at the LHC. The conclusion 
about the dominance of the hard diffractive mechanism over the e.m. one, when light nuclei 
are used, also holds in the RHIC kinematics. 

Our results are also valid for the diffraction in resolved photon-nucleus interactions. Since 
in this case several other effects are also important, we will discuss hard diffraction in 7 A 
interactions in a separate publication. 

II. SUPPRESSION FACTOR FOR HARD PROTON-NUCLEUS DIFFRACTION 

The derivation of the expression for the amplitude of hard diffraction in hadron-nucleus 
scattering combines features of soft coherent diffraction in hadron-nucleus scattering and 
hard coherent diffraction in DIS on nuclear targets. Therefore, we shall briefly review co- 
herent soft and hard diffraction below. 



A. Soft coherent proton-nucleus diffraction 



At high energies, the cross section of soft coherent hadron-nucleus diffraction (diffraction 
dissociation), a D A Dl can be economically and reliably calculated using the Glauber-Gribov 



multiple scattering formalism 
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(color) fluctuations in the projectile 
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2 ll| g eneralized to include the so-called cross section 
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hA 
DD 



j d2h i^j daP(a)\T A (b,a)\ 2 -^J daP{a)T A (b,a) 



(1) 



In Eq. ([I]), b is the impact parameter (the two-dimensional vector connecting the trajectory 
of the projectile with the center of the target nucleus); P(cr) is the probability to find 
in the projectile a hadronic configuration, which interacts with target nucleons with the 
cross section a; T A is the projectile-nucleus scattering amplitude in the impact parameter 
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space [20] 

T A (b,a) = l-exp(~<rT A (b)) , (2) 
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where T A (b) = / dzp A {b,z) and p A is the nucleoli density normalized to the number of 
the nucleons A. The energy-dependence of a^ D is determined by the energy-dependence of 
P(cr), which is implied [161 ] . Note that in Eq. (j2j), we neglected the slope of the elementary 
hadron-nucleon scattering amplitude compared to the nuclear size and we assumed that the 
elementary scattering amplitude is purely imaginary, which is a good approximation at the 
high energies that we consider in this work. 

The function P(cr) describes the probability that the incoming hadron interacts with 
target nucleons with a given cross section a. In other words, -P(c) describes cross sec- 
tion fluctuations in the energetic projectile. The notion of P(cr) is introduced in order to 
have a compact phenomenological description of soft coherent diffraction in hadron-nucleon 
and hadron- nucleus scattering. As follows from Eq. (JT]), ignoring cross section fluctuation, 
i.e. setting P(cr) oc 5(a — cr tot ), would result in the unacceptable result that a 1 ^ = 0. 

The formalism of cross section fluctuations is based on the simple picture of diffractive 
dissociation in the laboratory reference frame developed by Feinberg and Pomeranchuk 22] 



and by Good and Walker [23{]. In this picture, the incoming hadron is represented by a 
coherent superposition of eigenstates of the scattering operator. Since different eigenstates 
correspond to different a, the scattered state is in general different from the incoming state, 
but it has the same quantum numbers. This corresponds to the process of diffractive disso- 
ciation. 

One should note that the formalism of cross section fluctuations implicitly uses the as- 
sumption of the completness of the scattering eigenstates and, hence, it is applicable only 
at t ~ 0. At t ^ 0, the diffractive final state can be produced as a result of some effec- 
tive interaction or as a result of hard parton scattering (for sufficiently large t), which have 
nothing to do with the cross section fluctuations in the projectile. 

The function P(cr) is different for different projectiles (protons, pions, photons). For the 
proton, P(a) has a narrow dispersion around a = a tot , where cr to t is the total proton- nucleon 
cross section. Therefore, one can Taylor-expand the integrand in Eq. ([1]) around a tot and 
keep only first two non- vanishing terms 

a hA « J d 2 bT 2 {b) e -a tot T A (b) (3) 



In Eq. ([3]), LO a is the energy-dependent parameter, which is proportional to the proton- 
proton diffractive dissociation cross section and which controls the magnitude of cross section 
fluctuations 

uj„ = 4 Ms - 1 . (4) 



[f daaP(a 

Equation ([T]) can be interpreted as follows. The incoming proton diffractively dissociates 
on the front face of the target nucleus. The corresponding scattering amplitude squared 
is proportional to u a of ot T\(b). On the way through the nucleus, the produced diffractive 
state interacts with all nucleons of the target and becomes partially absorbed (suppressed). 
The corresponding soft suppression factor can be read off from Eq. (j3J), 

T soft = ex P (-°tot t a(P)) . (5) 

Note that since we have assumed that the dispersion of P(&) around a = er tot is small, the 
soft suppression factor depends only on cr tot . 



B. Hard coherent diffraction in DIS on nuclear targets 



Inclusive and coherent DIS on nuclear targets measure usual and diffractive nuclear parton 
distribution functions (PDFs), respectively. The theory of leading twist nuclear shadowing 
of usual and diffractive nuclear PDFs is based on the Gribov's relation between nuclear 
shadowing and diffraction 
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S 
13,18 



24|, Collins' factorization theorem for hard diffraction in 
4 ] and the QCD analysis of the HERA data on hard diffraction in DIS on hydrogen j^, 
9j, see 18( for the review and references. 
According to this approach, the nuclear shadowing correction, 5fj/A, to the nuclear PDF 
of the flavor j, fj/A = Afj/ N — Sfj/A, is expressed in terms of the proton diffractive PDF 

,D(4) 
Jj/N 
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24] 



Sxf j/A (x,Ql) 



(1 - ir) f 



1 + T] 2 

x PA (b, z x ) PA {b, z 2 ) e fc*"w(*i-«) 



oo rxo 
2 



(i 2 b I dz x j dz 2 J dx F f3f^\x,Ql,xp,t = 0) 



(6) 



In Eq. (E]), x and Qq are the Bjorken variables; xjp is the target longitudinal momentum 
fraction loss or the longitudinal momentum fraction carried by the diffractive exchange 
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(the "Pomeron"); (3 = x/xjp; x = 0.1; a J eS is the effective rescattering cross section of 
the intermediate diffractive state; 77 is the ratio of the real to the imaginary parts of the 
elementary diffractive amplitude. In the present analysis, 77 = tt/2 (op(0) — 1) = 0.185 Q, 9}. 

It is important to mention that since the t-dependence of the nuclear form factor is much 
steeper than that of the nucleon diffractive structure function, it is a good approximation to 
use f^\x,Ql,xjp,t min ) [t min « -(xpttin) 2 } instead of f?$\x,Ql,x F ,t) in Eq. ©. More- 
over, in the considered range of Bjorken x, t m \ n is small enough so that f^w (x, Qq, xjp, t m in) 
and f^/N (x,Ql,Xp,t = 0) practically coinside. 

As follows from Eq. (jSJ), a 3 eS determines the nuclear correction to Sfj/A due to the in- 
teraction with two and more nucleons (the interaction associated with the rescattering of 
the intermediate diffractive state). This cross section is defined in terms of the nucleon 
diffractive {ffffl and usual (f j/N ) PDFs \m,m\ 

< S (^Ql)= 7 , n2 J dxM.j3fffi(p,Cfi,x F ,t = 0). (7) 
xfj/ N (x, y J J x 

Figure [T] presents a J eS for gluons and w-quarks as a function of x at On = 4 GeV 2 . We used 
the recent QCD analysis of HI data on hard diffraction at HERA Q, \$ and CTEQ5M fit to 
inclusive PDFs 25]. 



Also, since Eq. (jTJ) involves the nucleon diffractive PDFs at t — 0, one has to make 



an assumption about the t-dependence of the nucleon diffractive PDFs. Experimentally, 
the t-dependence of the diffractive cross section is found to be practically constant as a 
function of /3, while the contribution of the gluon diffractive PDF increases strongly with a 
decrease of/? js), {J. Hence, in this work, we assume that all PDFs have the same exponential 
t-dependence, 



f^\{3,Qlxjp,t) = e-^f^'^Qlx^t = 0) , (8) 

where B^g = 6 GeV -2 is taken from the recent HI measurement with the leading proton 
spectrometer jsj. 

Figure [2] presents an example of our calculation of nuclear shadowing for nuclear PDFs of 
40 Ca and 208 Pb as a function of Q 2 and Bjorken x. The solid curve corresponds to the ratio 
fj/A/ (Afj/N) for u-quarks; the dotted curve corresponds to gluons. The lower set of curves 
corresponds to Q 2 = Ql = 4 GeV 2 . In addition to nuclear shadowing given by Eq. (jBj), 
we have introduced an enhancement (antishadowing) of nuclear gluon PDF on the interval 
0.03 < x < 0.2, which is modelled by requiring the conservation of the momentum sum rule, 



-B diB \t\ f D ^)i 



see e.g. 



18j. The two other sets of predictions for Q 2 = 10 GeV 2 and Q 2 = 100 GeV 2 are 



obtained by NLO DGLAP evolution. 

We would like to point out that the numerical analysis of nuclear shadowing presented in 



this work differs from our earlier analyses, see e.g. 



because we now use the most recent 



HI fits to nuclear diffractive PDFs and a different value of B^g, see Eq. (JHJ). However, 
the changes in the predicted nuclear shadowing are not large. The ratio fj /a/ {Afj /n) f° r 
w-quarks in Fig. [2] is very similar to our earlier result 18] . In the gluon channel, the ratio 

18]. 



fj/A/(Afjw) i n Fig. [2] is similar to the lower-gluon-shadowing scenario of 

Next we turn to nuclear diffractive PDF. In the Glauber-Gribov approach, the nuclear 
diffractive PDF of^the flavor j, is expressed in terms of the proton diffractive PDF 



ffl^ as follows 



xf°f{x,Qlxjp) = 4 7r/5/^ 4) (x,Q 2 ,x JP ,t = 0) J d 2 b 



dz e 



lQ7if^\x,Qlxjp,t = 0) I d 2 b 



' r r'" v z e -< r eff( x >Qo)/ 2 J!T dz 'pA(b,z') p A ^ ; z) 

_ e -<r> eS (x,Q 2 )/2T A {b) 



Ij/N 



(9) 



The last line is an approximation valid at small xjp, when the effect of the coherence length 
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FIG. 2: The ratio of the nuclear to nucleon PDFs, fj/A/(Af j/N ), for 40 Ca and 208 Pb at Q 2 = 4, 
10 and 100 GeV 2 as a function of Bjorken x, see Eq. (|6|). The solid curve corresponds to u-quarks; 
the dotted curve corresponds to gluons. 



[the factor exp(ixjpm N z)} can be neglected. In the opposite limit of large x F , xp > 0.05, the 
dominant contribution to the nuclear diffractive structure function is given by the impulse 
approximation, i.e. by Eq. Qj where a 3 eS is set to zero. In Eq. ([9]), the superscripts (3) and 
(4) denote the dependence on three and four variables, respectively. Note that similarly to 
Eq. ([3D, we neglected the slope and the real part of the elementary diffractive amplitude in 
Eq. ©. 

One can quantify nuclear diffractive PDFs by introducing the probability of diffraction 



for a given parton flavor j, Pi 



diff 



24) 



diff " xf^Ql) ■ m 

An example of the evaluation of the probability of hard diffraction in DIS according to 
Eq. (TTU1) is presented in Fig. [31 where Pj iff is plotted at fixed Ql = 4 GeV 2 as a function of 
Bjorken x. The left panel corresponds to DIS on 40 Ca; the right panel corresponds to DIS 
on 208 Pb. For comparison, the results for DIS on hydrogen are also given by thin curves. 
The solid curves correspond to the M-quark channel; the dotted curves correspond to the 
gluon channel. 
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FIG. 3: The probability of hard diffraction in DIS, P%:a, see Eq. (jlOp . The left panel corresponds 
to DIS on °Ca; the right panel corresponds to DIS on 208 Pb. For comparison, the results for DIS 
on hydrogen are given by thin curves. The solid curves correspond to the u-quark channel; the 
dotted curves correspond to the gluon channel. 



We would like to point out the following two features of P^ m presented in Fig. [3j First, 
the difference between the quark and the gluon channels is very small. While the quark 
and gluon diffractive and usual nuclear PDFs are different, their difference cancels to a large 
extent in the ratio P^ m (the cancellation is larger for heavier nuclei). Second, even for such a 
heavy nucleus as 208 Pb, P^ iS < 0.36, which should be compared to the asymptotic (A — > oo 
and <7g ff — > oo) upper limit Pj ifT = 0.5. An examination shows that while close to the center 
of the nucleus (b « 0), the probability of diffraction is very close to 1/2, the contribution of 
the nuclear edge significantly dilutes -Pjiff- 

Equation (Q can be interpreted as follows. The incoming virtual photon fluctuates into its 
hard diffractive component long time before the photon interacts with the target (we ignore 
the effect of the finite coherent length). The hard diffractive component elastically rescatters 
on the target nucleus, which gives the suppression factor [1 — exp(— ^a J eS (x, Qq) T^(6)] 2 , and 
emerges as the final hard diffractive state. 

One should note that the approximate expression for [the last line in Eq. (j5J)] 
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FIG. 4: Feynman graphs representing the first two terms of the multiple scattering series ([9]) for 
the nuclear diffractive parton distribution ff/^ ■ 

corresponds to the first term in Eq. ([3]) since the elastic contribution to DIS is absent 
(suppressed by the smallness of a e . m ). Therefore, the analogy between Eqs. (j3J) and (J9]) 
enables us to introduce the attenuation factor characterizing the suppression of hard coherent 
diffraction in DIS on nuclear targets due to nuclear shadowing, 

T& = e*j>{-o> B (x i C&)T A (b)) . (11) 

The graphical representation of Eq. ([9]), when only the interaction with one and two 
nucleons of the target is retained, is shown in Fig. HI The right graph helps to understand 
why T^ ai ^ is driven by a J eS . As seen from the graph, the strength of the hard rescattering is 
defined by the XN — > XN cross section summed over all X. This cross section is nothing 
but the ratio of the 7*iV — > XN to the 7*iV — > X cross sections summed over X, which 
corresponds exactly to a J eS defined by Eq. (J7J). 

It is important to emphasize that, in general, the calculation of T h ^ rd is model-independent 
only for the interaction with one or two nucleons. For the interaction with N > 3 nucleons, 
we implicitly used the so-called quasi-eikonal approximation in Eq. (fTTT) , which assumes that 
the diffractively produced state elastically rescatters on the nucleons. This approximation 
is equivalent to the observation of the small dispersion of -P(cr) used in the derivation of 
Eq. ©. 
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C. Hard coherent proton- nucleus diffraction 



As an example of hard coherent diffractive processes on heavy nuclear targets, we consider 
the hard coherent diffractive production of two jets in the reaction p+A — > 2 jets + X + A In 
this process, A denotes the nucleus; X denotes the soft diffractive component; the invariant 
mass of the jets provides the hard scale. 

The cross section of the p + A — > 2 jets + X + A reaction can be readily obtained by 
generalizing the well-known expression for the dijet inclusive cross section in hadron-hadron 
scattering [26] and by introducing the new quantity, the screened nuclear diffractive PDFs 

?D(3) 
Jj/A > 

( ^3 (T p+yl^2jets+X+ J 4 

dxi dp\ dX]p 

oc /</p(^i.^® 8 H^.^.^)D^(v^*oi 2 i^-, (12) 

id, _ 

k,l=q,q,9 

where fi/ p are the usual proton PDFs; — > kl)\ 2 are the invariant matrix elements 



for two-to-two parton scattering given in Table 7.1 of 26]; x\ and X2 are the light-cone 
momentum fractions of the proton and the nucleus active quarks, respectively; px is the 
transverse momentum of each of the final jets; Q e g is the effective hard scale of the process. 
For the simplification of our analysis, we consider only the case of 90° hard parton scattering 
in the center of mass, which constrains function of x% = fixjp) and Ql s 

*i = -^, <& = 4p&, (13) 
pxjp s 

where \/s is the proton-nucleon invariant energy. The term "screened PDF" means that 
this parton distribution contains certain soft suppression effects, i.e. the screened PDF is 
suppressed compared to the analogous PDF extracted from hard processes. 

The derivation of the expression for the screened nuclear diffractive PDFs, f^^i ls carried 
out similarly to the derivation of Eq. (jHJ) [see also Fig. [5] 

ff / f(x,Ql,x 1P )^Arrf^\x,Ql,x 1P ,t = 0) J <PbT${b) e -(«^CH-ir(-.^W , (14) 

where ffff is the screened diffractive PDF of the nucleon, which enters the QCD description 
of the p + p — > 2 jets + X + p reaction; <7 to t is the total proton-nucleon cross section; a 3 eS is 
the effective rescattering cross section of Eq. ([7]). In Eq. ffT4l) . we neglected the slope and 
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FIG. 5: Feynman graphs representing the first two terms of the multiple scattering series for the 
p + A — > 2 jets + X + A scattering amplitude. 

the real part of the elementary p + N — > 2 jets + X + N scattering amplitude and a small 
longitudinal momentum transfer in the p + N —>■ 2 jets + X + N vertex. 

It is important to emphasize that in the case of hard coherent proton-nucleus diffraction, 
the nuclear suppression factor, T h ^ d , is a product of the soft and hard suppression factors 
introduced previously, 

rppA _rppArp~f*A flr\ 

-'hard ~~ -'soft -'hard • \ LO ) 

This can be understood from Fig. [51 which represents the single and double scattering 
contributions to the p + A — > 2 jets + X + A scattering amplitude. The rescattering cross 
section of the middle graph is cr to t; the rescattering cross section of the right graph is a J eS 
(we assume that all diffractive intermediate states correspond to the same rescattering cross 
section). Therefore, the resulting nuclear attenuation, which results from the sum of the 
middle and right graphs, is driven by the cr tot + cr J eS cross section. 

Equation (|T4j) can be interpreted in two complimentary ways. On the one hand, one 
can start from soft diffractive dissociation of protons on heavy nuclei, see Eq. (J3j) . Since we 
are interested in the hard diffractive component of the diffractive dissociation cross section, 
one has to take into account the additional suppression of nuclear diffractive PDFs given 
by T^ av ^. As a result, one arrives at Eq. ffT5l) . On the other hand, one can start from 
the expression for inclusive diffraction of protons on nuclei, which is proportional to the 
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FIG. 6: The suppression factor A J of Eq. (|16p as a function of Bjorken x at Q 2 = 4 GeV 2 in the 
LHC (y/s = 8.8 TeV) and RHIC (\/s = 200 GeV) kinematics. The solid curves correspond to the 
u-quark; the dotted curves correspond to the gluons. 



nuclear diffractive PDFs Since the final diffractive state contains a soft component, 
which is partially absorbed by the nucleus, one should take into account this suppression by 
introducing the factor Tj^, which represents the probability of the absence of soft inelastic 
interactions at a given impact parameter b. 

We quantify the suppression of the nuclear screened diffractive PDFs compared to 
the nucleon screened diffractive PDFs J^yjy by introducing the factor A J 

f j/N ■{PiQiBr) fj/ N (P,Q 2 ,xp) J 

= B m (4ir) J d 2 bTl(b) exp [-(a tot (s) + ^ ff (x, Q 2 ))T A (b)] , (16) 

where in the last line we introduced the slope of the t-dependence of the screened diffractive 
PDFs, 

fffi(P,Cf,z F ,t) = exp(-B diS \t\) f$ ) (P,Cf,z F ,t = 0) . (17) 

Note also that f^\/3, Q 2 , x P ) = J dtf°ff((3,Q 2 ,XiP,t). 

Certain features of Eq. ([IS]) deserve a discussion. First, while the diffractive PDFs depend 
separately on (3 and xp, the suppression factor A J depends only on their product x = fixjp 
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FIG. 7: The suppression factor A J of Eq. (|16|) as a function of the atomic number at x = 10 3 and 
at Q 2 = 4 GeV 2 in the LHC (vs « 9 TeV) and RHIC (s/s = 200 GeV) kinematics. The labeling 
of the curves is the same as in Fig. [6J 



in our approach. Second, at the LHC energies, where a to t is of the order of 100 mb, the 
dependence of A J on a 3 eS is rather weak. Therefore, we expect that A- 7 is rather similar for 
different parton flavors j. In addition, since the slope B<he is independent on the hard scale 
Q 2 in our approach, A- 7 has very weak dependence on Q 2 , which enters only through the 
Q 2 -dependence of a J eS . 

In our numerical analysis of Eq. ffTBT) . we used the following input. The slope of the t- 
dependence of was taken from the recent HI measurement of hard inclusive diffraction 



- a. 



in DIS on hydrogen, B diS = 6 GeV 



The total proton-nucleon scattering cross section, er tot , was taken from 27J] 



a tot (s) = 21.7 s - 0808 + 56.08 S - a4525 . (18) 
The effective cross section crL was evaluated using Eq. (HI and the recent QCD fits to 

nn 

the HI measurement of hard inclusive diffraction on hydrogen [8, 9], see Fig. [TJ 
For the nucleon density pa, the two-parameter Fermi model was used 28]. 
Figures [6] and [7] present the results of our calculations of A- 7 in the LHC and RHIC 

kinematics. The LHC kinematics corresponds to \/s ~ 9 TeV per nucleon for proton- 
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p 



FIG. 8: The ultraperipheral p + A — > X + A scattering. 



nucleus collisions 29( ; the RHIC kinematics corresponds to y's = 200 GeV. Figure [6] gives 
A J as a function of Bjorken x at Q 2 = 4 GeV 2 . The solid curves correspond to the w-quark; 
the dotted curves correspond to the gluons. Despite the fact that A- 2 is of the order of 
unity at the LHC and of the order of several units at the RHIC energies, the corresponding 
suppression of hard diffraction is very large because in the absence of the suppression, 
nuclear diffractive PDFs are enhanced compared to the nucleon diffractive PDFs by the 
factor ff®lff$ oc A*. 

Figure [7] presents the ^4-dependence of A- 7 at x — 10 -3 and at Q 2 = 4 GeV 2 , i.e. at fixed 
<7g ff . As seen from Fig. [7J the A-dependence of A J is rather slow. A simple fit gives that 
A J oc A - 28 at the LHC and RHIC. 



III. HARD DIFFRACTION AND ULTRAPERIPHERAL PROTON-NUCLEUS 
COLLISIONS 

In proton-heavy nucleus ( 208 Pb, for example) collisions, most of the diffractive events 
(~ 80%) will be generated by the scattering of the proton off the coherent nuclear Coulomb 
field at large impact parameters, p + A — > p + + A X + A [lij, see Fig. These 
ultraperipheral proton-nucleus collisions open a possibility for studies of hard photon-proton 
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interactions at extremely high energies and allow one to probe the gluon density in the proton 
at the values of Bjorken x, which are a factor of ten smaller (for the same virtuality) than 



those probed at HERA [30|, M, [32]. 



In this Section, we estimate the ratio of the jet production in hard coherent proton- 
heavy nucleus diffraction (p + A — > 2 jets + X + A) to the production of hard jets by 
the photon-proton interaction, where the photon is coherently produced by the elastically 
recoiled nucleus, p + A — > p + '-f + A — ■> 2 jets + X + A This corresponds to the situation when 
the generic final state X in Fig. [8] contains a hard two-jet component and a soft remaining 
part X. 

Qualitatively, we expect that the ratio of the hard dijet production due to these two 
mechanisms will be rather small because of the following two suppression effects. 

First, hard diffractive dijet production is suppressed by the factor discussed in Sect. [Til 
Second, the shapes of the parton distribution in the photon and in the screened nuclear 
diffractive PDFs are rather different. In the photon case, the dominant contribution to the 
photon PDFs comes from (3 ~ 1 corresponding to the kinematics where a pair of jets is 
at the rapidities close to the gap. On the other hand, in the screened nuclear diffractive 
PDFs at large virtualities, which are relevant for the measurements at the LHC, the main 



contribution comes from small (3, see Fig. 3 of [17]. 

Using the definition of the suppression factor A J (fT6l) . the hard coherent proton- nucleus 
diffractive dijet cross section can be written as [see Eq. ( [121 ] 

^3 a p+A^2jcts+X+A 

dx\ dp\ dx jp 
ocr h M*i,Q«)X(PxM>,C&)f?^ 

*7T Jl — 1 + Okl 

k,l=q,q,g 

In the second line of Eq. ( fl9|) . we introduced an additional suppression factor r^, 

f^=r if $\ (20) 

which, according to the discussion in Sect. [TJ takes into account the significant factorization 
breaking in hard hadron-hadron diffraction. 

In our numerical analysis, we used the following model for the suppression factor 

75 / f os C° \ _1 

r h = r^r = 0.75 / dxjp / dtf P/p {x P ,t)) . (21) 

iV l s J \Jl.5/s J-oo J 
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This expression is based on the phenomenological model of [lCj, which describes the sup- 
pression of diffraction at the Tevatron (y/s = 546 and 1800 GeV) by rescaling the Pomeron 
flux, fp/ p (xp, t), by the factor N(s). In Eq. ( f2~Tl) . the Pomeron flux is given by the following 
expression 

fr/pixp, t) = x l+2e+2a't "jf^ ' ( 22 ) 



where e = 0.1; a 1 = 0.25 GeV -2 ; (3ppp{t) is the Ppp form factor [101 ]. 

We also introduced the additional factor 0.75 in Eq. fT2Tl) in order to phenomenologically 
take into account the observation that the effects of factorization breaking should be larger 
in the elementary diffractive PDFs at t = [see Eq. ( Tl6|) ] than in the t-integrated diffractive 
PDFs [see Eq. (12T])]. 

The application of Eq. (|2ip at the RHIC and LHC energies gives 

r h = ^L, RHIC, 

rh = i^o ' LHC - (23) 



Note that as follows from the definition (1211 . the suppression factor rh is assumed to be 
x jp-independent . 

Next we discuss the hard coherent dijet production in proton-nucleus scattering via the 
e.m. mechanism, when the nucleus coherently emits a quasi-real photon which interacts with 
the proton and diffractively produces two hard jets, p + A-^p + j + A — > 2 jets + X + A, 
see Fig. [HJ The corresponding cross section can be written as a sum of the resolved and 
direct photon contributions (the separation into the resolved and direct components is only 
meaningful in the leading-order calculation) 

k,l=q,q,g 

+ £ h/M,Ql s ) 1 ^5{P-l)Y J \M{;H^kl^-^-, (24) 

' Xjp — 1 + Okl 

i,3, _ 
k,l=q,q,g 

where n(xjp) is the flux of equivalent photons 30| expressed in terms of xp instead of the 
photon energy u (note the factor 1/xp coming from the 1/um the spectrum of the equivalent 
photons); fj/ y is the PDF of the real photon; X] 1-^(^7 y ^(llf are invariant matrix elements 
for the direct photon-parton scattering, see Table 7.9 in 26|; r em . is a phenomenological 
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factor describing the factorization breaking for the resolved (hadron-like) component of the 
real photon. The exact value of r e . m . is uncertain: It ranges from r e . m . = 0.34 [33] to r c m . « 1 
with large errors [34J . Since our analysis is a simple leading-order estimate, we conservatively 
take r e m . = 0.5. 

The flux of equivalent photons approximately equals 30j 

n(xjp) » In ( — -1 ) , (25) 

7T \RAXpPUhJ 

where Z is the nuclear charge; 7 is the Lorentz factor (7 « 3000 for pP6 scattering at 



the LHC 



291]); i?^ = 1-145 A 1 / 3 fm is the effective nuclear radius; pi a b is the momentum 



of the nucleus in the laboratory frame (pi a b ~ 2.75 TeV for pPb scattering at the LHC). 
In practice, we used a more precise formula for the flux of the equivalent photons, which 
reduces the result of Eq. f[2"5"j) by 11% [3^]. 

We are now ready to estimate the ratio of the hard diffractive dijet cross sections corre- 
sponding to the hard and e.m. mechanisms, R, 

J3^p+A^2jcts+X+A / j3 p+A^2jcts+X+A 

= dxMp I Z Mf • < 26 > 

where the involved cross section are given by Eqs. (fT2"|) and (1241) with the equal coefficients 
of proportionality. In the simplified kinematics that we use, at given p? and xjp, the ratio 
R depends only on (3. 

We considered two cases: The dijet production summed over gluon and quark jets and 
the production of two heavy-quark jets (c and b quarks). The resulting ratios R at pr = 5 
GeV and xjp = 10~ 4 , 10 -3 and 10 -2 as functions of j3 are presented in Fig. [91 The left panel 
corresponds to quark and gluon jets; the right panel corresponds to heavy-quark jets. 

The results presented in Fig. |9] deserve a detailed discussion. The dependence of the ratio 
R on xjp is not too strong and can be explained as follows. The main contribution to the 
xjp-dependence of R at fixed (3 comes from the changing of x\. As xjp is decreased, X\ is 
increased, which diminishes the role played by the gluons in the projectile. As explained 
in the following, it is the gluon contribution that increases R. Hence, R decreases with 
decreasing xjp. Note that the dependence of diffractive PDFs on xjp, f®$(f3,xjp,Ql s ) oc 
l/x 1 p 2e , see Eq. (1221) . is similar to the 1/xp ln(l/xjp)-behavior of the e.m. cross section. 
Therefore, these two factors weakly affect the x p- dependence of R. 
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FIG. 9: The suppression of hard diffractive dijet (quark and gluon jets) production compared 
to e.m. coherent dijet production in proton-Pb scattering at the LHC. The suppression factor R 
of Eq. (|26p at pr = 5 GeV and xjp = 10 -4 , 10 -3 and 10 -2 as a function of (3. The left panel 
corresponds to quark and gluon jets; the right panel corresponds to heavy-quark jets. 



The dependence of R on (3 is much faster and reflects different shapes of the proton 
diffractive PDFs and PDFs of the real photon. While the proton diffractive PDFs times 
(3 are flat in the (3 — » limit, the photon PDFs times (3 grow. This explains why R 
approaches zero when (3 is small. In the opposite limit, (3 — > 1, diffractive PDFs vanish and 
the e.m. contribution wins over due to the non-vanishing direct photon contribution: R — ► 
as /3 — > 1. 

In Fig. [HI the ratio R at its peak is much larger for the production of quark and gluon 
jets than for the production of heavy-quark jets. An examination shows that this effect is 
due to the large gluon diffractive PDF, which in tandem with the large gg — > gg hard parton 
invariant matrix element 26], works to increase R in the presence of the gluon jets. 



We used the following input in our numerical analysis of the ratio R. We used the LO 
parameterization of the real photon PDFs from Ref . 36| . We have also checked that the use 
of a different parameterization 37| leads to rather similar predictions. 

For the nucleon diffractive PDFs, we used the recent QCD fit to the HI data on hard 
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FIG. 10: The suppression factor R of Eq. (|26p for proton- 40 Ca scattering at the LHC. The labeling 
of the curves is the same as in Fig. [9] 

inclusive diffraction in DIS on hydrogen j^, 9]. The suppression factor A J , which enters 
Eq. (|26|) at the scale Q 2 = Q 2 ea = Ap 2 ^ = 100 GeV 2 , was evaluated using Eq. (|T6l) with 
a 3 eS (x,Q 2 ) at the same Q 2 = Q 2 cS = 100 GeV 2 scale, see Eq. 

The (5-function for the direct photon contribution was numerically modeled in the follow- 
ing simple form 

6{p-l) = l^-±p-^, with e = 0.01. (27) 

It is instructive to examine how our predictions for the suppression factor R change, 
when the heavy nucleus of 208 Pb is replaced by a lighter nucleus of 40 Ca. Note that for pCa 
scattering at the LHC, y^i = 9.9 TeV and 7 ~ 3700 [29]. We expect that the ratio R will 
significantly increase because of the reduction of the flux of the equivalent photons [the flux 
is proportional Z 2 (125}) ]. 

Figure [10] presents the ratio R for 40 Ca. The labeling of the curves is the same as in 
Fig. [9] As can be seen from the comparison of Figs. [10] and [9] the ratio R increases by the 
factor « 7 — 10, when going from 208 Pb to 40 Ca. 

Besides the LHC, RHIC also has a potential to measure hard diffraction in proton-nucleus 
scattering. We consider a typical example of the corresponding RHIC kinematics with 250 
GeV protons scattering on 100 GeV/per nucleon nuclei (the corresponding y/s ~ 320 GeV 
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FIG. 11: The suppression factor R of Eq. (|26p in the RHIC kinematics and at px = 5 GeV and 
xp = 5 x 10~ 3 as a function of (3. The solid curves correspond to quark and gluon jets; the dotted 
curves correspond to heavy-quark jets. 

and the Lorentz dilation factor is 7 ~ 100). Producing sufficiently high diffractive masses, 
e.g. M\ = 500 GeV 2 , one accesses the typical kinematics of hard diffraction, xp = 5 x 10 -3 
and (3 > 0.3. Note also that the suppression of hard diffraction at RHIC is approximately 
four times smaller than at the LHC, see Eq. fl23l) . 

We studied the suppression factor R of Eq. (I26I) in the considered RHIC kinematics at 
Pt = 5 GeV. The resulting values of R as a function of j3 are presented in Fig. [TTJ The solid 
curves correspond to quark and gluon jets; the dotted curves correspond to heavy-quark 
jets. 

As seen from Fig. [TTJ, the factor R at RHIC is larger than at the LHC. This is mostly a 
consequence of the decrease of the flux of equivalent photons when going from the LHC to 
the RHIC kinematics. 

IV. CONCLUSIONS 

Using the Glauber-Gribov multiple scattering formalism and the leading twist theory of 
nuclear shadowing, we developed a method for the calculation of coherent hard diffraction 
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processes off nuclei. We showed that soft multiple rescatterings lead to the factorization 
breaking of hard diffraction in proton-nucleus scattering, which is larger than the well-known 
factorization breaking of diffraction in hadron-nucleon scattering. 

Based on these results, we compare the hard diffractive to e.m. mechanisms of hard 
coherent production of two jets in proton-nucleus scattering. We study the xjp, ft and A- 
dependence of the ratio of the dijet production cross sections due to the two effects, R, at 
the LHC and RHIC kinematics. We separately study the case when the final jets consist of 
quarks and gluons and the case when the final jets consist of heavy (c and b) quarks. 

Our results can be summarized as follows. For proton- 208 Pb scattering at the LHC, hard 
diffraction is suppressed compared to the e.m. contribution, especially at xp = 10~ 4 and 
large ft, e.g. ft > 0.8, see Fig. The suppression is very strong for the production of heavy- 
quark jets, see the right panel of Fig. [91 The physical reason of the suppression is the strong 
coherent Coulomb field of 208 Pb, which enhances the e.m. mechanism of hard diffraction. 

Replacing 208 Pb by 40 Ca, the hard diffractive mechanism becomes compatible to the 
e.m. one in the case of the production of quark and gluons jets, see the left panel of Fig. [TU1 
However, like in the case of 208 Pb, the production of heavy-quark jets is dominated by the 
e.m. mechanism, see the right panel of Fig. [TOj 

As a result of the smaller Lorentz dilation factor 7 at RHIC, the factor R at the RHIC 
kinematics is larger than at the LHC. 

Our results suggest the following experimental strategies. First, the use of heavy nuclei 
in p A scattering at the LHC will provide a clean method to study hard real photon-proton 
scattering at the energies exceeding the HERA energies by the factor of ten. Second, taking 
lighter nuclei and choosing the appropriate kinematics, where the e.m. contribution can be 
controlled, one can effectively study the factorization breaking in nuclear diffractive PDFs. 
Third, in the same kinematics, a comparison of the dijet diffractive production to the heavy- 
quark-jet diffractive production will measure the nuclear screened diffractive gluon PDF. It 
can be compared to the nuclear diffractive PDFs, which will be measured in nucleus-nucleus 
ultraperipheral collisions at the LHC. 
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